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Abstract- Modelling and manufacturing of radar absorbing 
material are proposed using multilayer composite nanostructured 
materials and the recently introduced winning particle 
optimization algorithm. The study concerns tile of materials with 
dimensions of 0.3mx0.3m made of nanostructured composite 
materials, which consist on epoxy-resin and industrial grade 
carbon nanotube (CNTs) fillers. The industrial grade CNTs were 
appositely chosen for their low costs, in order to be applied in 
great amount to build large tile of composite materials. Here 
modelling takes into account for an extended frequency band (5 to 
18 GHz), for several incidence angles of the electromagnetic field 
(0 to 80°), and for the minimization of the electromagnetic 
reflection coefficient. At last simulations are compared with 
measurements of reflection coefficient. Despite some errors mostly 
due to the manufacturing process, simulations are in good 
agreement with measurements, showing an interesting approach 
to design multilayer radar absorbing materials. 

Keywords- Radar Absorbing Materials; Evolutionary 
Computation; Winning Particle Optimization; Modeling; Carbon 
Nanotubes; Composite Nanostrucutred Materials; Layered 
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I. INTRODUCTION 

Radar absorbing materials (RAM) acquired importance 
since their first applications in military field [1-6]. Germany 
pioneered the first aircraft to use RAM during World War n, in 
the form of the "Horten Ho 229", in order to reduce the radar 
signature, and concerned with radar camouflage for submarines, 
developed "Wesch" material, a carbonyl iron powder loaded 
rubber sheet about 0.3 inches thick and a resonant frequency at 
3 GHz. Composite materials considered in this work are based 
on epoxy matrix reinforced with carbon nanomaterials. These 
latter were chosen taking into account their good microwave 
absorption behaviour [7] and the low market prices. In fact, the 
economic aspects, normally neglected in small laboratory 
applications, are important in real applications where the 
amount of carbon nanopowders fillers could be relatively high. 
In such scenario, a good compromise in terms of 
cost/performances was obtained using industrial grade 
multiwall carbon nanotubes (MWCNTs) with a cost of about 
300 $/kg. MWCNTs with commercial name of NANOCYLTM 
NC7000 (diameter around 9.5 nm, length 1.5 \xm, purity 90%) 
were bought at NANOCYL. Bi-component Epoxy-resin 
PrimeTM 20LV (density 1.123 g/cm 3 ) with Hardner (density 
0.936 g/cm") was used as matrix. RAM is based on layered 
materials where each layer is made of carbon nanostructured 
composites materials in different weight percentages. The 
recent developed winning particle optimization search 



algorithm (WPO) was applied in order to design and optimize 
multilayer materials able to effectively absorb microwaves in 
the range 5 GHz to 18 GHz. The paper is organized in 5 
distinct sections. In Section II, mathematical model of absorber 
is shown; in Section III, WPO is described; in Section IV, 
numerical design of RAM is discussed and in Section V, RAM 
manufacturing and electromagnetic performance comparisons 
between simulation and manufactured RAM are shown. In this 
section the adopted experimental setup is accurately described 
showing the NRL-arch test fixture. It is crucial to highlight the 
interdisciplinary research activity through nanomaterials, 
electromagnetic wave propagation theory, composite materials 
manufacturing techniques and search computation algorithms. 
All of them are required to design the "quasi perfect absorber". 

II. RAM MATHEMATICAL MODEL 

WPO algorithm [8-9], has been applied to solve a quite 
complex problem consisting in the design and optimization of 
microwave absorbing multilayer materials. 

Optimized RAM is based on nanomaterials in particular 
carbon nanostructured composites materials. Conductive fillers 
have been uniformly dispersed in an epoxy resin at different 
weight percentages (0.5, 1, 2, 2.5, 3 wt %). In Table I, six 
different composite materials and the index code adopted by 
WPO to identify each material are shown. 



table I 

CODING OF MATERIALS IN THE DATA BASE 



Materials 


Materials Code 


Epoxy- Resin 


1 


MWCNT,0.5 wt% 


2 


MWCNT, 1.0 wt% 




MWCNT,2.0wt% 


4 


MWCNT, 2.5 wt% 


5 


MWCNT, 3.0 wt% 


6 



The resulting composites samples have been dielectrically 
characterized in a previous work [10-11] by waveguide 
measurements to recover the dielectric properties of the 
composite materials in the database (DB). In this paper the DB 
of these dielectric parameters as function of frequency is used 
to supply the WPO modeling algorithm. In Table n, a brief 
summary of the dielectric permittivity of composite materials 
at 6GHz, 12 GHz, 18 GHz, is shown. After that, numerical 
simulations of RAM have been carried out. In particular, 
electromagnetic analysis has been performed integrating the 
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forward/backward propagation formalism to the in-house built 
WPO, thus able to carry out optimization upon oblique 
incidence angles over a finite angular range. 

TABLE II DIELECTRIC PERMITTIVITY OFCOMPOSITE NANO STRUCTURED 
MATERIALS 



Materials 


Real Part and Imaginary Part 


Fl= 6 Ghz 


F2= 12 Ghz 


F3=18 Ghz 


Epoxy- Resin 


3.11 


0.10 


3.07 


0.08 


3.07 


0.08 


MWCNT,0.5Wt% 


5.52; 


0.90 


5.33 


1.81 


5.21 


2.83 


MWCNT, 1.0Wt% 


9.52 


5.27 


6.75 


2.88 


5.39 


1.95 


MWCNT, 2.0 Wt% 


11.42 


6.89 


8.82 


4.48 


7.45 


3.38 


MWCNT,2.5Wt% 


12.90 


11.28 


10.86 


10.35 


9.71 


9.79 


MWCNT, 3.0Wt% 


18.21 


14.50 


17.29 


13.65 


16.72 


13.31 



The developed code minimizes reflection and transmission 
coefficients under the thickness minimization constraint. 
Finally, broadband quasi-perfect absorbers in the band 5-18 
GHz are achieved combining the filler families, i.e., exhibiting 
a loss factor (LF%) greater than 90% in most of the band, for 
thicknesses ranging between 5 and 12 mm. Main goal of design 
and optimization work was to achieve values (< -10 dB) of 
reflection coefficients (RC), for angular apertures within 40° in 
the most part of the frequency band. The design/optimization is 
basically a minimization procedure which seeks the best trade- 
off between structure thickness (to be minimized) and absorbed 
EM power (to be maximized). The absorbing power has been 
quantified through the so-called "loss factor" LF%, defined as 
follows [10-11]. 



LF% = (] 



l-\RC' 



\TC\ 



! )-100. 



(1) 



In RAM systems the transmission coefficient TC=0 
(because of the perfect electric conductor (PEC) layer at the 
end of the multilayer), and the RC is the WPO-optimized 
reflection coefficient. 

Both are expressed in the linear form. LF% physical 
meaning is related to the fraction of the incident power 
vanishing inside the materials because of localization and 
dissipation phenomena. Since the calculation of RC becomes a 
crucial point for RAM optimization goal, a rigorous approach 
for computing reflected waves at each layer has been adopted 
[11-13]. 

Both normal and oblique incidences have been evaluated 
during the automatic optimization. In Fig. 1, a general scheme 
of the electromagnetic absorbing multilayer structure is 
presented. Each layer is generally denoted by indexx 
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Fig. 1 General Multilayer scheme of electromagnetic absorbing structure; in 
RAM sy stemsthe back layer is a perfect electric conductor (PEC). 



As far as geometric properties of multilayer structures are 
concerned, thickness of each layer can range from to 10 mm, 
while the number of layers is upper-bounded to maximum 10 
layers. The effective number of layers is established by the 
WPO optimization procedure which is able to remove one or 
more useless layers from the multilayer structure. The order of 
layer materials is not a priori fixed; instead WPO algorithm 
decides which material need to be used for each layer. 

In RAM applications, RC at the air-absorber interface can 
be evaluated by the following equation relating the free space 
impedance to the input impedance seen at the air- multilayer 
structure interface 



RC' 



ry TM ITE _ 7 TM ITE 
^TM ITE + ^TM/TE 



(3) 



RC™ /TE =20Lo gi 



ry TM I TE _yTM ITE 
ry TM I TE ,ryTM ITE 



where Zo = 377 O is the free space impedance and Zi is the 
input impedance at the first air-absorber interface. The input 
wave impedance of the multilayer, backed by PEC, can be 
expressed iterating the following equation for x-th-layer [11- 
13]. 



ryTM ITE _ TM ITE 



Z™/-cos(Mj + M™™sin 



n™ ,TE cos(kj^ + jzZ ,TE Mk x t x ) 



(4) 



In (4), t x = y x _! -y x is the x-th layer thickness expressed in 
m, whereas the wave number k x of each layer along y axes is 

kl= (2tt/) 2 Mo«oM a e„ cos 2 3 X = (5) 

= faff M £ MrM* - j £ ri) COs2 9 X 

where f is frequency of the incident electromagnetic wave 
in Hz,s' is the real part of x-th-layer permittivity and e" is the 
imaginary part of layer permittivity. 

Complex waves appear in oblique incidence and lossy 
dielectrics problems [13]. Because of the wave number become 

complex-valued, e.g., k = ft- ja , the angle of refraction and 
possibly the angle of incidence may become complex-valued 
too. In calculating k x by taking square root of (5), it is required, 
in complex-wave problems, to get the correct signs of their 
imaginary parts, such that evanescent waves are described 
correctly. This leads to define an "evanescent" square root as 

follows. Let e rx = [s rx - js r ^ with s rx >0 for an absorbing 

medium, then 



kl = 



Ufa fYMo e o." a (e a - 
l -j^-( 2 ^f) 2 M £ M r . 



js r Jcos & x M s rx >0 (6) 
fc)cos 2 ^,if s rx =0 



Since all nanostructured composite materials available in 
the DB are with s n > then k x is as usually computed. 

HI. WINNING PARTICLE OPTIMIZATION 

WPO is a simple algorithm where at each time epoch of 
evolution, the particle, which best fit the objective function, is 
deputed to pilot the trajectory of the remaining particles within 
the multidimensional space of solutions, i.e., variables to be 
optimized. In Fig. 2, WPO flow chart is shown. At the 
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beginning of the iterations, particles are randomly distributed 
within the n-dimensional search space. 



Start WPO Algorithm 



Initial random particles 
generation 



Particles P, fitness evaluation 



Mark the best particle {BP") with 
(/-index ranging from 1 to max 

particles number 
I 



Extraction of all m-coordinates 
of each particles in the 
n-dimensional space 




l:n 






+ PK 







Mutation operation for some 
randomly position coordinates m 
of the particles p k ^ pg'i 



Check of boundary conditions 
for each single m-coordinate 




End WPO Algorithm 



Fig. 2 WPO flow chart 

At each iteration time, the objective function of each 
particle is computed. Then, particles are compared to each 
other using value of their objective function and the best fitting 
particle is marked with its proper index. Each particle position 
is completely defined by its coordinates and their number 
represents the dimension of the space where all particles jump 
searching for the optimal solution. Trajectory of each particle, 
except that of the best fitting particle, is defined in the 
following way: if mth-coordinate of a certain particle P k is 
lesser then the corresponding mth-coordinate of the current best 
particle PB q , then the new mth-coordinate will result from the 
sum of current particle and best particle mth -coordinates; while 
if mth-coordinate of a certain particle P k is greater than the 
corresponding mth-coordinate of the current best particle PB g , 
then the new mth-coordinate will result from the difference 



between the current particle mth-coordinates and the best 
particle mth-coordinates. Here, q is the index of the best 
particle found in each iteration and the condition which needs 
to be grant is q^k, i.e., all kth particles can be displaced except 
the best qth particle. Equation (7) describes the mentioned 
approach. Particle position at the ( i+1 jth iteration as a function 
of particle position in the previous (th iteration is. 



P*(i + l)=SR m P*(i)+ gPB« m (i) 



(7) 



where 

g=+l if P*(i)<PB« m (i) 

g=-\ if P*(i)>PBl(i) 

k = Particle index in the range [1 , Particle Number] 

q = Best Particle index in the range [1 , Particle Number] 

and k q 
m = ranges in [1 , n] 
n = Space Dimensions 
R m = random number in the interval [0,1] 
8 = Convergence Parameter 
i = current iteration in the range [ 1 , N] 
N = Number of iteration 

R m is used to randomly define the amount of displacement 
of each single coordinate. Such randomization develops the 
searching ability by conferring some causality to the amplitude 
ofjumps. Convergence parameters helps the final convergence 
of the WPO algorithm as follows. At the beginning of the 
iterations, the WPO must be capable to explore the largest 
space of solutions as possible. On the other hand, close to the 
end of the iterations, high amplitude of jumps could represent a 
drawback since a position representing the optimal solution 
could forbid. In order to improve the WPO convergence, a 
mechanism providing a progressive reduction of "maximum 
jumps amplitude" has been introduced and it is identified by 
the 5 parameter. Calling N the total number of iterations and 
supposing ith the current iteration, then the value of 5 is given 
by 



5 = a 1 



N + l 
j(R m -0.5) 



(8) 



(9) 



where J is the "jump amplitude" defined by the user as a 
numerical parameter in the main program of WPO, in this 
scenario it has been valued to 3 after several trials to assure a 
good performance in domain solution exploration. 

In (8), the meaning of a and S parameters is in fast or slow 
convergence of the WPO algorithm In particular, the higher 
the a, the greater the initial jumps amplitude and the greaterthe 
distances where particles will start to explore the search space. 
As far as S is concerned the higher its value, the faster the 
convergence to a suboptimal solution. The parameter a is made 
using J and /? m such that (R m - 0.5) is a random number 
between [-1/2, +1/2]. The randompart of a develops thesearch 
capability of the algorithm. Even though mutation operation 
has not been applied on the subsequent reported optimizations 
(mutation=0%), it has to be noticed that mutation can be 
included in WPO algorithm to further increase the random 
search capability. At the end of each WPO iteration, checking 
of boundary conditions for the new set of particles coordinates 



IJMSCI Vol.2 lss.2 2012 PP. 31-38 www. ij-ms.org©World Academic Publishing 



~ 33 ~ 



International I ournal of Material Science 



IIMSCI 



is computed in order to avoid overcoming the constraints . The 
entire loop cycle can be iterated several times in order to get 
the required objective function minimization. 

IV. RAM DESIGN AND OPTIMIZATION 

The presented modeling procedure, takes into account for 
several incidence angles and, the design/optimization 
procedure can be forced to run for a one particular incidence 
angle or for a range of angular incidences where the global 
objective function has to be optimized. In particular the 
incidence angles here range between 0° and 80° with step 10°. 
As far as layer materials are concerned each one can be chosen 
among all of different manufactured composite materials. In 
Table I, all the composite materials available in the DB are 
reported and each composite material has been coded using an 
integer number. 

WPO can access to DB in order to allocate the most 
appropriate materials for each layer of the multilayer RAM. In 
WPO procedures the modeling functions are called by the main 
program to compute: materials intrinsic impedances, in- 
outgoing refraction angle, TE/TM layer wave impedance, 
reflection coefficient (RC), and the loss factor (LF). The 
multilayer absorbing structures considered in the design 
procedure can have up to 10 layers so dimensionality of layer 
space is 10 i.e., m - 10, 9. . . 1. 

Since the optimization are for both: layer thickness and 
layer material type then WPO algorithms need to be structured 
in order to take into account for both quantities. 

In WPO algorithm we have: T k (i) = [**(/) t k (i), , t\ (/)]. 

as array of thicknesses, where T k {i) is the thickness of m-th 

layer, and P k (i)=[p k w (i\p k (i\ ,p k (i)] > as array of material 

type, where (/) is the type of material associated to the mth 
layer. 

The equations related to layer thicknesses and layer 
material for each single particle are simply updated according 
to, 

T k (i + l) = S R ml T k (i) + g mckness TBI (0 ( 10) 

P k (i + 1) = round(s R m2 P k (i) + g Material PB« m (ij) 

Where, R„a and R„q are random numbers in the closed 

range 0-1, T k (i) is the current ith-particle, mth-layer- 

thickness at the rth-iteration, while P k (i) is the current kth- 

particle, mth-layer-material at the (th-iteration. 

Rounded values in material equation are required since the 
index of material of each layer must be an integer number 
ranging from 1 to the number of materials available in the DB 
of materials. 



(11) 



+i if J i(0<raj(i) 

"1 if r*(i)>7Bj(/) 

W W =+1 if Pl(})<PB q m (i) 

g Material = 'I $ *£(0 >«£('") 



where q is the index of best particles TB and PB; k is the 
indexof particles with k^q; n is the dimension number; m is the 
current dimension ranging between [l,n]. In the presented 



implementation, the value of convergence parameter S is set to 
1. 

WPO try to minimize a global objective function (GOF) 
which in turn is built using elementary objective functions 
(EOF). The EOF for TM and TE modes are named as 
CostRC TM , CostRC je, for reflection coefficient. The formal 
definition of the EOF for TM and TE modes are shown in 
(12),(13),(14). We can observe that for each particle (Pa), the 
corresponding EOF is evaluated over the entire frequency band 
and over the entire incidence angular range, adopting 
frequency and angular steps chosen by the user before starting 
the WPO, 

freq=fmaf Q=9 max 

£ ^RCjPaJreqj) 

freq-f min \$-#mm 



CostRC TM (Pa) 



(12) 



CostRC TE {Pa)-- 



freq=fmsa/ $=& max 

£ ^RCjPaJreq^) 



freq-f min \ #=#min 



where Pa is the current particle, freq is the frequency step, finin 
and finax is the frequency band start and stop, is the current 
angular step, 9 min and 6^ represent the angular range bounds. 
The definition of the EOF for thickness is 

CostT{Pa) =[t w (Pa) + t 9 (Pa) + .... + t l (Pa)] (13) 

where t m is the thickness of the m-th layer. 

A weighting factor called a ranging in (0<a<l), weight 
CostT w.r.t. CostRC. Such weighting factor is chosen by the 
user and its meaning have to be intended as the capability of 
the tool to design the multilayer structure making privilege to 
the electromagnetic performances w.r.t. the thickness when gis 
close to 1. Final GOF is a linear combination of the described 
EOFs, 

GOF(Pa) = a OFl(Pa) +(!-«) OFl{Pa) ( I 4 ) 



where 



OF\(Pa) = 
OF2(Pa)~- 



, (Pa) CostRC T1 



(Pa) 



( CostRC TU y 
'I A 
~ CostT(Pa) ~ 
B 

A and B are normalization factors: A is the product 
between the frequency step number and the angle step number, 
whereas B is the product between the maximum layer number 
and the maximum layer thickness. In the following subsection 
two RAM models are shown. 

A. RAM First Simulation 

In the first RAM simulation a weighting factor 
a=0.9 returned the scheme of layered material shown in Fig. 3. 

The electromagnetic wave impinges upon the first layer 
made of epoxy -resin. Incidence and reflection angle is called 0°, 
ranging from 0° to 80°. Total thickness of multilayer RAM is 
11.77 mm and the final number of layers is 4. Tile is 
composed on the top by the lowest lossy material (epoxy -resin) 
as the first layer up to the highest lossy material in the DB 
(MWCNTs 3 wt%) as third layer backed by PEC. This 
structure can be viewed as a graded lossy system in the first 
three layers and a resonant structure in the fourth layer, where 
the chosen materials are not the highest lossy materials 
(MWCNTs 0.5% wt). In particular graded lossy systems are 
composed by a quasi-wave impedance matched material where 
the first interface between free space and epoxy-resin is the 
best coupling material able to assure the lowest reflection 
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coefficient at the first air-absorber interface. MWCNTs lwt% 
and MWCNTs 3wt% instead have lower characteristic wave 
impedances and greater losses with respect to epoxy-resin and 
as a consequence they are placed in sequence in the layers 
order. 

Incidence Wave \ l 4 Reflected Wave 



Infinite length 




-11.77 mm 

Fig. 3 RAM multilayer simulated: Resin -epoxy, MWCNTs lwt%, 
MWCNTs 3wt%, MWCNTs 0.5 wt%,PEC 

In Fig. 4, reflection coefficients TE/TM (dB) are shown. It 
can be noticed that for normal incidence angles, TE and TM 
curves coincides whereas around 1,9, 12, 15 GHz, there are 
some incidence angles where curves show values of reflection 
coefficient TM even lower with respect to that at normal 
incidence. Such angles are denoted as Zenneck angles for lossy 
materials. 
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Fig. 4 RAM Reflection Coefficient (dB) TE andTM 
The electromagnetic absorbing capability of a certain 
material is mainly due to two simultaneous conditions: 
"dielectric losses" within the composite material and 



"impedance matching condition" which in turns is a function of 
thickness of the absorbing layer. The impedance matching 
condition takes place when the electromagnetic wave 
impedance at the first interface between air and RAM is close 
to that of free space i.e., about 377 Q. This is why in the 
simulation the first layer is logically made of epoxy-resin 
which possesses the lowest electric permittivity and so the 
higher electromagnetic wave impedance. The electromagnetic 
wave impedance at the first interface is computed using (4) and 
depends on dielectric parameter of all materials in the 
multilayer. The dielectric losses are described by the imaginary 
part of the permittivity (s") while, thickness of the material is 
connected to the microwave wavelength within the material. 
The higher the dielectric losses, the greater are the 
electromagnetic power dissipation phenomena in the RAM 
composite material. Unfortunately, usually, the higher the 
dielectric permittivity, the lower is the electromagnetic wave 
impedance, which in turn mismatch the impedance matching 
condition with free space [11]. This is why, the RAM is 
composed of different layers able to progressively match the 
wave impedance while simultaneously assuring the right 
dissipation and absorbing behaviour. In Fig. 5, Loss Factor 
TE/TM is shown. Loss Factor (%) remains above 90% for 
incidence angles up to 40°. Such values mean that RAM is able 
to absorb the most part of the electromagnetic field reflecting 
back only a very little amount of the incoming energy. The 
extended frequency band 6 to 18 GHz where loss factor is 
greater than 95%, shows the interesting theoretical 
performance of this RAM. 




10 12 14 

Frequency (GHz) 




10 12 14 
Frequency (GHz) 

Fig. 5 RAM Loss Factor (%) TE andTM 
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In Fig. 6, the asymptotic convergence of WPO evolution is 
shown. Plot of GOF is shown. The simulation is made on 800 
iterations and 256 searching particles. The time required to 
complete the computation strongly depends on the amount of 
particles and iterations and by the computer performances in 
terms of random memory allocation and processing speed. 
After 600 iterations the behaviour of curve is almost 
asymptotic showing a probably sufficient convergence of the 
WPO algorithm in searching the best solution to the RAM 
design and optimization problem. The computing time also 
depends on the way in which code is written. Here Matlab has 
been used. Some procedures were parallelized but others 
require nested loops of cycle to take into account for particles, 
frequencies, and incidence angles. These are very time 
consuming. 




200 300 400 500 600 700 
Iterations of WPO Algorithm 

Fig. 6 WPO evolution and convergence 



B. 



RAM Second Simulation 
In the second RAM simulation, a lower weighting factor 
a=0.6 returned the scheme of layered materials shown in Fig. 7. 
The multilayer RAM is composed of three different layers, the 
first one, where the electromagnetic wave arrives, is epoxy- 
resin; the second one is made of MWCNTs at lwt%, the third 
layer is made of MWCNTs at 2wt%. The total thickness is 
about 5.72 mm. In Figs. 8 and 9 the electromagnetic 
performances of RAM are shown. 



Incidence Wave 



Reflected Wave 



Infinite 




= 2.7958 mm 
= 0.4132 mm 
:2.5176 mm 



~5.72 mm 

Fig. 7 RAM multilayer simulated: Resin-epoxy MWCNTs 1 wt%, 
MWCNTs 2wt%, PEC 
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Fig. 8 RAM Reflection Coefficient (dB) TE andTM 
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Fig. 9 RAM Loss Factor % TE and TM 
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The graphs show two peaks of absorption at 0° incidence 
angle for TE/TM around 7 and 16 GHz. When the incidence 
angle grows to 50°, 60° and 70° in the TM incidence mode, the 
first absorption packs shift to 8, 10 and 11 GHz while the 
second peak disappear. 

V . RAM MANUFACTURING AND ELECTROMAGNETIC 
PERFORMANCE MEASUREMENT 

The second simulated RAM multilayer scheme has been 
chosen as the model to try of building a large RAM tile. Each 
layer was manufactured over the previous using the sample 
holder as shown in Fig. 10. The manufacturing procedure is 
available in [9]. In Fig. 11, the manufactured RAM tile is 
shown. 




Fig. 10 Sample holder to build multilayer RAM 




respect as much as possible the RAM scheme and thickness 
performed by WPO simulation, in the case of Fig. 7. 



Incidence Wave 



Reflected Wave 



30 cm 




~6.1 mm 

Fig. 12 RAM multilayer manufactured: Resin-epoxy MWCNTs 1 
wt%, MWCNTs 2wt%, PEC 
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Fig. 13 RAM multilayer manufactured: resin-epoxy MWCNTs 1 wt%, 
MWCNTs2wt%,PEC 




Fie. 11 Final manufactured RAM multilavertile 
The practical realization of the large RAM tile was not easy 
since the dimensions 30x30 cm required a lot of nanomaterials 
and the related dispersion process in the epoxy-resin resulted 
hard. Another very difficult task was to respect the thickness of 
each layer because the manually homemade deposition of each 
layer above the other one caused errors in the value of 
thickness and in the homogeneity of layers along the surface. 
The great manufacturing problems rose with the layer made of 
MWCNTs at 2wt%. In fact, such composite doesn't remain 
enough viscous so tile manufacturing process is quite difficult, 
especially for respect thickness constraints. At the end, the 
final manufactured multilayer RAM is lightly different with 
respect to the initial simulation reported in Fig. 7. In Fig. 12 the 
manufactured RAM multilayer tile scheme is shown. The total 
thickness of tile is about 6.1 mm. The simulation of such new 
tile is shown in Fig. 13. Observing the simulation it can be 
highlighted that even small variations of thickness layers are 
able to affect the reflection coefficient lowering the absorption 
properties of RAM. This is why it is very important to try to 



To evaluate the electromagnetic reflectivity, a NRL arch 
measurement (bistatic reflection method), has been used [14]. 
In this configuration, two antennas are placed for transmitting 
and receiving signals respectively and the microwave 
reflectivity at different incident angles can be measured. Using 
this configuration, the reflection properties of materials, can be 
measured for different incident angles. It should be noted that 
in bistatic reflection measurements, the reflection is dependent 
on the polarization of the incident wave. Incident waves with 
parallel and perpendicular polarization usually result in 
different reflection coefficients. Besides, special calibration is 
needed for free-space bistatic reflection measurements [14]. In 
Fig. 14 bistatic measurement system and the tile in composite 
material is shown. Bistatic measurement system is here based 
on Agilent software 857 IE (material measurement), and 
Agilent PNA-L vector network analyser. Antenna used are 
SATIMO Dual Ridge Horn SH2000 covering the frequency 
range 2-32 GHz. After the calibration of the NRL bistatic 
system, the reference in reflection coefficient has been taken 
adopting only metal plate i.e., without RAM tile. After that 
measurement of a known sample consisting in ECCOSORB 
AN73 absorber [7], has been performed in order to be aware 
about errors in the NRL measurement setup. Confidence of 
measurements was within 2 dB of interval with respect to 
reflection properties declared in the ECCOSORB data sheet. 
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Fig. 14 Bistaticmeasurement system located at (www.saslab.eu) 

Measurement of the manufactured RAM tile using 
incidence angle of 0° has been performed for TE and TM mode. 

In Fig. 14 comparison between simulation of the reflection 
coefficient of the new manufactured RAM tile and real 
measurements is shown. The angle of electromagnetic wave 
incidence is 0°. It can be noticed that the simulated curve 
shows similar behaviour with respect to the measured one. 

Errors and differences in simulated and measured reflection 
coefficient are mainly due to: 

• difficulties in manufacturing RAM multilayer tile 
where each layer thickness ranges of about 0.5mm between its 
minimum to its maximum value; so the reported value 
represent an average value. To be aware of the manufacturing 
errors, in Fig. 14 we show the measurements of RAM called 
TE_0, TM_0, TE_90, TM_90. The numbers and 90 mean 
that the RAM tile placed on the metal plate was rotated 
respectively of 0° and 90°. It can be noticed that there are some 
differences in the measured reflection coefficients which can 
be explained by errors in the manufacturing process. 

• accuracy in determination of permittivity of composite 
materials performed using wave guide method. 

• accuracy in free space measurement of reflection 
coefficient using bistatic method. 

V. CONCLUSION 
In this work modelling, optimization and manufacturing of 
radar absorbing materials have been discussed. The design and 
optimization took place using the recently introduced winning 
particle optimization search algorithm which has been applied 
to optimize the electromagnetic absorbing properties while 
contemporarily minimizing the overall thickness of multilayer 
material. Following the simulation, a large tile of radar 
absorbing material has been manufactured by using industrial 
grade of multiwall carbon nanotubes and measurement of 
electromagnetic reflection coefficient performed by NRL-arch. 
Comparison between simulation and measurement show 
interesting agreement. Nevertheless some errors occurred in 
manufacturing process which affected electromagnetic 
performances of manufactured radar absorbing material. 
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